The diffuse neutrino flux from supernovae: upper limit on the electron neutrino 
component from the non-observation of antineutrinos at SuperKamiokande 
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I derive an upper bound on the electron neutrino component of the diffuse supernova neutrino 
flux from the constraint on the antineutrino component at SuperKamiokande. The connection 
between antineutrino and neutrino channels is due to the similarity of the muon and tau neutrino 
and antineutrino fluxes produced in a supernova, and to the conversion of these species into electron 
neutrinos and antineutrinos inside the star. The limit on the electron neutrino flux is 5.5 cm~^s~^ 
above 19.3 MeV of neutrino energy, and is stronger than the direct limit from Mont Blanc by three 
orders of magnitude. It represents the minimal sensitivity required at future direct searches, and 
is intriguingly close to the reach of the Sudbury Neutrino Observatory (SNO) and of the ICARUS 
experiment. The electron neutrino flux will have a lower bound if the electron antineutrino flux is 
measured. Indicatively, the first can be smaller than the second at most by a factor of 2-3 depending 
on the details of the neutrino spectra at production. 

PACS numbers: 97.60.Bw,14.60.Pq 



The observation of the diffuse flux of neutrinos from 
core collapse supernovae has become a concrete possibil- 
ity in the last years. If achieved, it would represent a 
milestone of neutrino astronomy, as it would allow faster 
progress with respect to searches for an individual super- 
nova signal, and would provide several tests of neutrino 
properties, of the physics of core collapse and of the cos- 
mological rate of supernovae. 

Currently, the best experimental result on the diffuse 
supernova neutrino flux (DSNz^F from here on) is the up- 
per limit on its electron antineutrino {9^) component, put 
by SuperKamiokande (SK) above 19.3 MeV of neutrino 
energy . This limit comes from the non-observation of 
events due to 9^ absorption on protons in water above 
the background. At 90% confidence level, it reads: 



'S>eiE > 19.3 MeV) < 1.2 cm 



(1) 



This bound approaches the range of theoretical predic- 
tions of the flux suggesting that a positive signal 
may be seen in the near future. Such signal will add 
to the information given by the data from SN1987A 

MM 

The effort to detect or constrain the electron neutrino 
{ve) component of the DSNi^F is at least as important as 
the search in the De channel. Indeed, it has the potential 
to give the very flrst data on neutrinos (instead than an- 
tineutrinos) from supernovae, with crucial first tests of 
theory - e.g., of neutrino transport in dense media - in 
the neutrino channel. In spite of its importance, how- 
ever, the search of the i^e component of the DSNi^F has 
had less progress than that of D^, . This is because in wa- 
ter an exclusive (one fiavor only) signal is not possible, 
and non-water experiments are forced to smaller volumes 
than SK, resulting in a lower sensitivity. The best exper- 
imental information on the v,, fiux is the bound: 



at 90% confidence level, found in 1992 at the Liquid 
Scintillator Detector (LSD) of Mont Blanc from the 
non-observation of absorption on ^'^C New, 
stronger bounds will come from the Sudbury Neutrino 
Observatory (SNO) O and from the upcoming experi- 
ment for Imaging Cosmic And Rare Underground Signals 
(ICARUS) M- 

In this paper I show that the i^e diffuse fiux is already 
strongly restricted by the SK result itself. Indeed, the 
limit on Pf, , Eq. , can be easily translated into an anal- 
ogous result for Vg. This new bound improves greatly on 
the LSD limit. Clearly, it has the same information con- 
tent as the direct bound (^3) , and so it is important mainly 
because it defines the minimum sensitivity required for 
future direct searches in the Ve channel. 

Briefly, what allows to place a limit on the i^e flux from 
the result on De is the combination of two elements. The 
first is the fact that the fluxes of muon and tau neutrinos 
and antineutrinos {i^^, Vt,v^,Vt) produced inside the star 
are equal to first approximation. The second element is 
that neutrinos undergo flavor conversion due to masses 
and mixings on their way from production to detection. 
This implies that at least part of the detected J/g ipe) flux 



originates from the produced and i>t 



Vr)- Their 



$e(25 MeV < £; < 50 MeV) < 6.8 ■ 10^ cm'^s" 



(2) 



origin from species that have similar fluxes at production 
provides the connection between the and channels. 
This same connection will also give a lower bound on 
the component of the DSNz/F if the Ve flux is detected, 
thus giving important guidance to direct searches of the 
i^e flux. 

Some generalities are in order. In a core collapse super- 
nova neutrinos are present as a thermal gas, right outside 
the core, from where they diffuse out carrying away about 
99% of the 0(10^'^) ergs of gravitational energy that is 
liberated in the collapse. This energy is shared among 
neutrinos and antineutrinos of the three flavors, e,^,r, 
in comparable amounts. 
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Before flavor conversion, the energy spectrum of each 
flavor of neutrinos is expected to be nearly thermal with 
average energy of 10-20 MeV The non-electron 

species, i'^, Vt, and v^., interact with matter only 
via neutral current exchange, which is flavor independent 
and, at the lowest order in momentum transfer, parity 
preserving [l^ . This implies that and Vt [v^ and 
Vr) have equal fluxes and so they can be considered as a 
single species, (i^x)- The spectra of fx and Dx are sim- 
ilar, with differences due to weak magnetism corrections 
to the neutrino-nucleus cross section. These cause the 
to have a ~ 7% lower average energy One also 

expects softer spectrum of Ve with respect to Pe-, because 
the interaction of (9^) with matter is dominated by 
capture on neutrons (protons) and matter is highly ncu- 
tronized deep inside the star. In summary, the hierarchy 
of average energies E'os Eq^ ^ £'oe ^ Eq^ is expected 
at the point of decoupling of neutrinos from matter, and 
is supported by numerical calculations (e.g. [la llSl 1211 ) ■ 
Indicative values of the average energies are: _Eoe ~ 9 — 16 
MeV, Ene - 12 - 18 MeV, E^x - 15 - 22 MeV. 

Between the production and the detection points, the 
neutrino fluxes are modified by two effects: the redshift 
of energy and the flavor conversion (oscillations) in the 
star and in the Earth. The latter is due to the interplay of 
neutrino masses, flavor mixing and coherent interaction 
of neutrinos with the medium (refraction). In cases where 
the conversion has no zenith nor energy dependence, the 
diffuse fluxes of v,, and of in a detector in a given 
energy bin Ei — E2 have the form: 



(1 -p)<i)x 



-f (l-p)0s , (3) 



where p (p) is the Ve (i^e) survival probability, i.e. the 
probability that a neutrino produced as i/g (Dg) is de- 
tected as Ve (f'e)- Here 0„, is the diffuse flux of i^^ 
(w = e,e,x,x) in the detector in absence of conversion. 
It depends on the redshift z, on the flux originally 
produced in a individual star, dNy^^E') / dE' , and on the 
comoving rate of supernovac, Rsn- 
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dN.^{E') RsN{z)dz 



dE' 



_ (4) 

(see e.g. [Tj). Here E and £" = (1 + z)E are the neutrino 
energy at Earth and at the production point; f2„j ~ 0.3 
and U,\ ~ 0.7 are the fraction of the cosmic energy den- 
sity in matter and dark energy respectively; c is the speed 
of light and Hq is the Hubble constant. 

Generally, the conversion of neutrinos depends on E 
and i?', on the matter density proflle along the neutrino 
trajectory - and thus on the arrival direction, for neu- 
trinos crossing the Earth - and on the neutrino masses 
and mixings. In particular, p and p are functions of the 
mixing angles 612 and ^13 (assuming the standard param- 
eterization of the mixing matrix, e.g. |2ll |'l and of the sign 



of the atmospheric mass squared splitting, Am 
9i2 angle is measured to be 
,2 



2 

31- 



tan^ e 



12 



0.3-0.73 



The 



(5) 



at 99.73% C.L. |2a|, while 6*13 is still unknown, with the 
bound sin^6'i3 <, 0.02 For brevity, I wiU refer to 

the intervals sin^ ^ 2 • 10~^ and sin^ ^13 ^ 3 • 10"''' as 
"small" and "large" ^13 respectively. The sign of Am§]^ 
is also unknown; the two possibilities are called "nor- 
mal hierarchy" (Am|j^ > 0) and "inverted hierarchy" 
(Am^j^ < 0) of the neutrino mass spectrum. 

Three scenarios of conversion in the star are most rel- 
evant here (see Table jilEilifij: 

(i) the "adiahatic normal", characterized by the nor- 
mal hierarchy and large ^13. In this scenario neutrinos 
undergo resonant conversion [2^ at matter density of 
~ 10"^ g • cm^"^. Their propagation is adiabatic, i.e. with- 
out transitions between the eigenstates of the Hamilto- 
nian, resulting in p = 0. The propagation of antineutri- 
nos is adiabatic too, with p = cos^ 6*12. 

(ii) the "adiahatic inverted" , with the inverted hierarchy 
and large ^13. Both neutrinos and antineutrinos prop- 
agate adiabatically, but the resonant conversion affects 
antineutrinos, so that one has p = sin^ 612 and p ~ Q. 

(iii) the "non- adiahatic" , characterized by small ^13. 
Here resonant conversion is made ineffective by the com- 
plete breaking of adiabaticity, and one has p = svc? 612 
and p = cos^ 0i2 for any hierarchy. 

In these three cases, p and p have very little dependence 
on the neutrino energy, thus justifying the use of the 
expressions For intermediate values of ^13 and nor- 
mal (inverted) mass hierarchy, things are as in case (iii) 
with the difference that p {p) is energy-dependent and 
< p < sin^ 6*12 (0 < p < cos^ 6*12)- The survival prob- 
abilities receive small corrections from oscillations in the 
Earth (Hill [111 and Shockwave effects The former 
will be discussed briefly later; the latter are neglected as 
they are smaller than 1-2% in the energy interval where 
the DSNi^F is largest, £: - 20 - 30 MeV 113 ■ I also 
neglect terms of order df-^ and higher in p and p. 

Now I derive the upper bound on the ly^ flux. $e- I con- 
sider high energy thresholds, Ei <; 20 MeV, motivated 
by the cuts necessary for background reduction at exist- 
ing i/g and i>e detectors [llll3l|. First, I flnd an expression 
that is always greater or equal than the ratio ^e/^s- 

$e ^ p0e + {l- p)(l)x 



(1 -p)(px 

This follows from Eq. Q by dropping the tpg term there. 
Next, I use the fact that (j)x/4>x < 1, because the i/^; spec- 
trum is slightly softer than the one, and has average 
energy at or below the threshold Ei (see e.g. fig. 4 in ref. 
01). I also consider that the ratio © is largest when p 
has its maximum value Pmax- Thus: 

''^-<^—\l+p(^-l)] . (7) 



1 - Pri 
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character 


liierarcliy 


P 


P 




adiabatic 


normal 





COS^ 6*12 


0.67 - 4.6 

( ^ A 
(< ^■") 


adiabatic 


inverted 


sin^ 612 





0.39 - 1 
(> 0.39 ) 


fully non-adiabatic 


any 


sin'' 9i2 


COS^ 012 


0.39 - 4.6 

( any ) 



TABLE I: The survival probabilities p and p, and the ratio 
"l?e/<l?e for different combinations of mass hierarchy and de- 
gree of adiabaticity of the neutrino propagation, under the 
conditions (j)e/4'x < 1 and <j)^/4'x < 1. The numbers in 
brackets give the results without those conditions. I used 
Ma,'x{(j}x / (fix) ~ 1.5. The small effects of the Earth matter 
and of Shockwaves, as well as ©(Sfs) terms are neglected. 

It appears that, in general, an upper limit on $e/$g de- 
pends on a variety of unknowns through the factor (p^/ (p^. 
However, the bound simplifies to 

1^ < -^4 , if p = or <t>J^r, < 1 . (8) 

From TablclHl have pmax = cos^ 6*12 + 0(10"^), with the 
second term representing the (positive) correction due 
oscillations in the Earth. From Eqs. (0) and ((SJ I 

get: 

|^^-47r + ^(10"')^4.6, (9) 
•Pg sni 6*12 

$e(-B > 19.3 MeV) ^ 5.5 cm^^s^^ , (10) 

if J? = or < 1. Explicitly, Eqs. hold 

under one or the other of the following conditions: 

(a) the adiabatic normal scenario is realized (large 6*13, 
p = 0). No other assumptions - on the neutrino fluxes 
or on the supernova rate, etc. - are necessary, and thus 
in this case the constraint H1U|I is robust. 

(b) any other oscillations scenario is realized, or the con- 
version pattern remains unknown, and we know that 
<j)el<t>x < 1- The latter condition will be called ^^Vx domi- 
nance" . It is generally valid at > 20 MeV for neutrino 
fluxes motivated by theory; I will return to this point in 
a moment. 

Eq. H10() was obtained using the smallest 6112 in Eq. (jSj) 
and includes the correction due to Earth effects. This 
was estimated numerically following refs. 0, for a 
large variety of neutrino spectra and supernova rate func- 
tions; it contributes to the final result (|1U|I by at most 
5%. About the significance of the bound ((Tn|l . using the 
likelihood information from refs. P, I find that the 
inequality <I>e/ sin^ 6'i2 < 5.5 cm~^s~^ holds at 98% con- 
fidence level. 

A stronger bound follows from Eq. (jSJ for the adiabatic 
inverted case {pmax — 0), with dominance: 

— < 1 , ^e{E > 19.3 MeV) < 1.2 cm^^g-i ^ i^n-j 
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FIG. 1: The ratio (j)e/4>x as a function of Eq^, for iSoe = 12 
MeV and ratio of and luminosities of 2,1,0.5 (uppper to 
lower curves). I parameterized the neutrino spectra and the 
supernova rate function as in |3, with = O-x ~ 2.3, a = 0, 
/3 = 3.4. The SuperKamiokande energy threshold E\ = 19.3 
MeV was used. 



The more general limit, Eq. H10|) . is the main result of 
this work. It is highly motivating for i^e detectors, as it 
approaches their sensitivities. These were estimated to 
be ~ 6 cm~^s~^ in the energy interval 22.5 — 32.5 MeV for 
SNO (using charged current ly^ scattering on deuterium) 
[l3| . and ~ 1.6 cm~'^s~^ in the interval 14 - 40 MeV for 
ICARUS (using absorption in liquid Argon) 0| . 

One may wonder how Eq. H10|l changes when one in- 
cludes the dependence of neutrino conversion on the en- 
ergy, on the redshift and on the arrival direction of the 
neutrinos. The derivation, and thus the results Q) and 
(jS)), still apply, provided that for p, p and pmax one uses 
their values averaged over energy, redshift and zenith an- 
gle, with the product of neutrino spectrum and supernova 
rate as weight, see Eq. Q. This procedure gives the cor- 
rection due to Earth effects, explaining the 0(10^^) term 
in Eq. 0. The limits JH)) and Ullfl remain valid, as they 
refer to limiting cases, where p is truly a constant. 

A second question is the validity of the condition of 
i/x dominance above the threshold Ei 20 MeV. While 
there are no direct tests of this, the condition is supported 
by current supernova simulations, which predict the I'e's 
to have comparable luminosity (within a factor of two) 
and a softer spectrum than and i/x- Fig. Ogives the 
ratio 4'e/4>x calculated using Eq. |0J with typical param- 
eters (see caption). One sees that the dominance is 
violated by at most ~ 20% when the component is 
very soft and suppressed in luminosity. 

What can be said on the Ve flux if evidence of the 
DSNi^F is not seen in the u,, channel but is found in 
the Vf, one? Similarly to what shown so far, a positive 
detection of can be translated into a lower limit on 
the Uf, component. Indeed, one can constrain the ratio 
<f>g/$e from above, using nearly the same steps as in 
Eas.l(^- (|ll|l . with the exchanges e <-> e, a; <-> i, p ^ p. 
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The analogous of Eqs. ® and 1)10(1 are: 

<; Max(0/0.) ^ 

$e C0S2 



1.73 Max 



n -1 



0(10" 



(12) 
(13) 



where the largest value of 612 from Eq. (jSJ was used. 
These results are valid for the adiabatic inverted case 
{p = 0) or for Vx dominance, (f>e/(t'x ^ 1- They depend 
on the maximum value allowed by theory of the ratio 
(jix/ipx, Max((/)j/0a;), just like the derivation of Eqs. © 
and (|10|l contains the maximum of (j>x/(t>x- The latter 
was set to 1 there. 

To estimate Max{(j)x / (f>x) requires a dedicated work, 
with the evaluation of the dependence of (f>x/4>x on the 
details of the Ux and Dx original spectra. For illustration, 
here I take Max{(f>x / (j>x) = 1.5, which is the largest ratio 
of the Px and Vx spectra found in rcf. [T^ in the energy 
interval = 20 - 80 MeV. With this value Eq. lO 
gives $e 0.39$e — O(10~^) (the minus sign indicates a 
negative Earth effect). 

For the adiabatic normal scenario {p = 0) with the Vx 
dominance one finds the analogous of Eq. Hll|l : 



< Max{(j)x/(l}x) 



(14) 



This gives $e k, 0.67$g for Max{(l)x / (px) = 1-5. 

The results of this paper are summarized in Table 
The Table helps the discussion that follows. 

What will we learn from a detection of the DSNi^F in 
the Ve channel, together with the limits l|10(l and ((11)1 ? 
If evidence is found above the value in (fTn|l . the physics 
and/or the data analyses involved would have to be re- 
vised, e.g. by relaxing possible priors or by reconsidering 
background subtraction. If both the Ve and the Pe data 
analyses are proved to be robust, one would have to ad- 
mit the possibility of a very large Vg component in the 
flux before neutrino conversion (i.e. 0e > (jix above the 
detection threshold), or invoke exotic physics that could 
give a larger i/g survival probability, Pmax (see Eq. ((TJ). 
A measurement above the limit ((11(1 would exclude the 
adiabatic inverted scenario with h'x dominance. A de- 
tection would also constrain the De flux indirectly. For 
example, Eq. 1(13(1 (with Max(02/0j,) = 1.5) tells that 
a measured i/g flux smaller than ^ 0.48 cm~^s~^ above 
19.3 MeV would be incompatible with a i/g flux as large 
as 1.2 cm^^s^^ in the same energy bin, thus improving 
the SK bound for 9^, Eq. (QJ. 

In summary, I found an upper limit of $e ^ 
5.5 cm^^s^^ above 19.3 MeV for the i^e component of the 
diffuse supernova neutrino flux, using the SK bound on 
the Ue flux, the properties of supernova neutrino spectra 
and neutrino oscillations. This new limit is the strongest 
available, improving the LSD bound by three orders of 
magnitude. It defines the necessary sensitivity for future 



searches of the Ve flux. Noticeably, the result proves that 
the Ve flux needs not be equal or smaller than the Vg one - 
as it would be the case without oscillations - but instead 
it can be larger by a factor of several above E ^ 20 MeV, 
depending on the oscillation scenario. It follows that, in 
spite of their smaller volumes, z^e telescopes can be com- 
petitive with SK in the search of the diffuse supernova 
neutrino flux. A dedicated analysis of the existing SNO 
data is compelling, in this perspective. 
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and M. Strassler for precious comments on the prelim- 
inary manuscript, and to M. Maltoni for useful discus- 
sions. I thank ICTP of Trieste for hospitality while work- 
ing on this paper, and acknowledge support from the 
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